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摘要 : 在 近 地 空间 目标 和 碎片 的 定 轨 预 报 任务 中 ， 需 要 准确 计算 热 层 中 性 大 气 密度 。 太 阳 辐 
射 约 27 d 周期 的 短期 震荡 是 影响 大 气 密度 的 一 种 重要 的 空间 环境 因素 ， 它 会 引起 大 气 密度 在 
全 球 尺 度 上 的 涨 落 。 选 择 CHAMP, GRACE-A 和 SWARM-C 这 3 颗 极 轨 卫 星星 载 加 速 仪 数 
据 推 导出 的 大 气 密度 资料 ， 提 取 其 中 的 27 d 短期 变化 信号 ， 与 太阳 极 紫外 辐射 (S10 指数 ) 的 
27 d 短期 变化 进行 多 元 回归 分 析 ， 研 究 热 层 大 气 密度 对 太阳 辐射 短期 变化 的 响应 规律 ， 及 其 
在 不 同 高 度 、 辐 射 水 平 、 昼 夜半 球 和 纬度 等 方面 的 差异 。 结 果 表 明 : 太阳 辐射 短期 变化 对 大 气 
密度 的 影响 与 高 度 负 相关 (主要 因素 )、 与 辐射 水 平 正 相 关 (次 要 因素 ); 对 白天 半球 的 影响 是 
夜间 半球 的 约 2 倍 ; 在 白天 半球 ， 辐 射影 响 随 着 纬度 增 大 而 减弱 ， 夜 间 半 球 恰 好 相反 。 将 以 
上 观测 结果 与 NRLMSISE00 模型 模拟 结果 对 比 ， 发 现 模型 低估 了 太阳 辐射 短期 变化 对 大 气 密 
度 影响 的 (50 负 ~60%)， 且 低估 了 高 度 差 异 ， 但 高 估 了 辐射 水 平 差异 。 利 用 SOHO 卫星 对 太阳 
26~34 nm 波长 极 紫外 辐射 1 h 分辨 率 的 测量 值 ， 研 究 了 大 气 密度 对 辐射 的 响应 延迟 时 间 ， 约 为 
18 h (0.75 d)。 研 究 结果 对 于 优化 大 气 模 型 建 模 、 修 正大 气 密度 具有 借鉴 价值 。 
关 键 词 : 空间 天 气 ， 热 层 大 气 密度 ， 太 阳极 紫外 辐射 ，27 d 短期 变化 
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1 引 B 


对 于 近 地 空 间 目标 和 碎片 的 定 轨 预报 任务 ， 热 层 中 性 大 气 密度 是 一 个 重要 的 输入 参数 ， 
会 随 着 日 地 空间 环境 、 轨 道 位 置 等 很 多 因素 变化 。 其 中 ， 太 阳 活 动 是 影响 大 气 密度 变化 的 最 
主要 能 量 源 ， 它 的 11 a 基本 周期 和 27 d 短 周期 震荡 都 驱动 着 大 气 密度 呈现 出 对 应 的 周期 变 
化 。Jacchia 和 Briggs” 较 早 地 发 现 1958 52 的 轨道 加 速度 存在 30 d 左右 的 周期 变化 ， 并 猜 
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想 这 与 太阳 辐射 的 27 d 周期 变化 有 关 ， 也 不 能 排除 与 地 磁 活 动 的 相关 性 ;随后 用 多 颗 卫星 
4 周期 变 率 证 实 了 27 d 太阳 辐射 对 大 气 密度 的 直接 影响 ”。 后 续 很 多 大 气 模型 都 按照 线性 
律 描述 大 气 密 度 与 太阳 辐射 的 关系 ， 并 从 测量 连续 、 工 程 可 用 性 角度 考虑 ， 将 射电 辐射 流 
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nt SS 


EJ. Fior 代 蔡 极 紫外 辐射 作为 模型 输入 ， 并 取得 较 好 的 应 用 效果 。 
人 们 也 关注 大 气 密度 对 太阳 辐射 变化 的 响应 时 间 ， 早 期 Romer ”发 现 大 气 密度 的 变化 


滞后 Fio.7 指数 约 1.06 d; Eastes 等 人 "发 现 中 午 的 大 气 密度 滞后 软 K 射线 变化 约 1.5 d; 


Guo 等 人 "分 析 CHAMP 的 资料 后 发 现 密度 相对 于 射电 辐射 存在 1 d 左右 的 延迟 。Ren 等 
人 “基于 物理 大 气 模型 的 模拟 发 现 大 气 密度 峰值 相对 于 极 紫 外 27 d 辐射 的 延迟 时 间 约 为 


0.9 d， 对 大 气温 度 的 模拟 显示 延迟 时 间 范 围 在 0.5~0.8 d o 


早期 研究 在 资料 空间 和 时 间 分 辨 率 上 存在 一 定局 限 ， 因 此 大 和 气 密度 短期 变化 规律 在 细 


节 上 研究 得 不 够 清晰 ， 例 如 ， 对 不 同 高 度 、 地 方 时 和 地 理 纬度 的 大 气 ， 辐 射影 响 有 无 差异 ? 
不 同 辐射 水 平时 ， 影 响 程 度 是 否 一 致 ? 为 了 回答 这 些 问题 ， 


本 文选 择 CHAMP、GRACE-A 


和 SWARM-C 等 3 颗 卫 星 的 高 分 辨 率 和 全 球 纬度 覆盖 大 气 密度 资料 作为 基础 ， 再 用 极 紫外 


辐射 指数 Sio 和 SOHO 卫星 26~34 nm 极 紫 外 波段 测量 值 直接 描述 太阳 辐射 的 27 d 变化 ， 


用 个 例 对 比 和 多 元 回归 分 析 等 方法 ， 研 究 大 气 密度 啊 应 在 不 同 高 度 、 纬 度 、 代 夜半 球 和 辐射 


水 平等 方面 的 细节 特征 和 延迟 时 间 。 这 些 研 究 对 于 优化 大 气 模型 建 横 ， 修 正大 气 密度 具有 借 
鉴 意义 ， 有 助 于 提高 近 地 空 间 目标 和 碎片 的 轨道 预报 精度 ， 也 能 为 高 层 大 气 变 化 物理 机 制 


完 提 供 观测 证 据 。 


2 数据 和 方法 


2.1 太阳 极 紫 外 辐射 的 测量 及 其 指数 


太阳 发 出 的 软 X 射线 、 极 紫外 辐射 (EUV) 在 穿越 地 球 热 层 大 气 时 被 吸收 ， 对 中 性 气体 
产生 加 热效应 ， 提 升 大 气 密度 和 温度 。 然 而 ， 正 因为 它们 被 热 
只 能 用 地 面 望远镜 测量 射电 波段 


测 ， 需 要 去 外 层 空间 测量 。 早 期 研究 缺乏 空间 卫星 手段 ， 


层 吸收 ， 导 致 在 地 面 上 无 法 观 


流量 近似 代替 极 紫外 辐射 水 平 “。1995 年 SOHO 卫星 发 射 升 空 ， 搭 载 的 “太阳 极 紫 外 监听 
器 ”(SEM) 开始 了 对 极 紫 外 辐射 的 常规 连续 观测 。SpaceWX 人 处理 了 波长 为 26~34 nm 的 测 
量 数据 并 归 算 成 与 射电 辐射 流量 指数 Fior 一 致 的 计量 单位 (sfu, 1 sfu=10-22 W-m-?-Hz~?), 


发 布 了 极 紫 外 辐射 指数 510。 该 指数 已 被 用 于 建立 Jacchia-Bowman 系列 模型 (JB2006 和 
JB2008)， 通 过 与 其 他 辐射 波段 的 指数 (Fio.7, Mio, Yio) 联合 使 用 研究 太阳 辐射 对 大 气 密 度 
影响 ， 取 得 了 很 好 的 效果 ””。 Sus 指数 的 分 辩 率 是 1 d， 对 于 工程 应 用 已 经 足够 。 然 而 ， 


为 研究 大 气 密度 对 辐射 变化 的 响应 时 间 ， 还 需要 更 高 分 辨 率 的 极 紫 外 辐射 测量 值 ， 本 文选 择 


SOHO 卫星 26~34 nm 波段 测量 值 1 h 分 辩 率 的 资料 。 
2.2 ”太阳 辐射 的 27 d 周期 震荡 
太阳 存在 25~35 d 的 自转 (赤道 地 区 自转 快 ， 两 极地 


X E 


转 慢 ， 即 较 差 自转 ” )， 因 此 


太阳 表面 的 若干 活动 区 时 而 朝向 地 球 ， 时 而 侧 向 ， 时 而 背 向 ， 引 起 到 达 地 球 的 辐射 强度 随 之 


3 期 EER, 


起 伏 ， 变 化 周期 与 太阳 自转 以 及 活动 区 所 在 纬度 
射 在 射电 波段 (Fio.7 指数 ) 的 27 d 震荡 周期 为 (27.18 +2.26)d 。 用 相同 的 方法 对 1997 一 
2022 年 的 极 紫 外 辐射 (Sio 指数 ) 进行 分 析 ， 统 计 得 到 震荡 周期 为 (26.85 2.11) d， 与 之 前 


的 结果 接近 。 
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TX. 我 们 曾 统 计 了 1960—2013 HAGA 


图 1 是 用 Sio 指数 


81 d 滑动 均值 )。 从 图 1 


27 d 震荡 最 显著 的 区 间 者 


处 于 太阳 活动 峰 年 附 


[14] 


绘 的 11 a 基本 周期 (81 d 滑动 均值 ) 和 27 d 短期 变化 (日 值 减 去 
a) 可 以 清晰 看 到 太阳 的 第 23、24、25 个 基本 周期 ， 图 1 b) 表明 
近 (2001—2005 年 ，2011 一 2016 F), ÆRE 


时 (2008 一 2009 年 ，2019 一 2020 年 ) 震荡 幅度 很 低 。 


T " T " T T T T TT 


— §,,(814) 
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图 1 ”资料 覆盖 时 间 范 围 内 的 太阳 11 a 长 期 变化 和 27 d 短期 震荡 
图 2 是 放大 局 部 时 间 后 的 SEM 测量 值 和 Sio 指数 ，a) 是 强 太阳 活动 期 2002 年 5 月 至 


2003 年 10 月 ， 


27 d 短期 变化 很 汪 


Sli; b) 是 弱 太 阳 活 动 期 2019 年 1 A 2020 4E 6 H, f 
然 也 能 看 到 若干 周期 变化 ， 但 极度 平静 时 (2019 年 7 月 至 12 A) 没有 显著 周期 变化 。 将 两 


段 时 间 资 料 进行 对 比 : a) 时 间 段 的 SEM 均值 是 b) 的 3.28 ff; SEM 标准 差 (反映 出 样本 周 
期 变化 的 离散 度 ) 在 a) 时 间 段 是 b) 的 12.5 倍 ， 可 见 太 阳 活 动 越 强 ，27 d 短期 震荡 越 显 著 。 


2.3 ” 热 层 大 气 密度 数据 


选择 CHAMP、GRACE-A 和 SWARM-C 共 3 颗 卫 星 的 加 速 仪 测量 资料 反 演 得 出 的 大 
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2 ”太阳 活动 峰 年 (2002 一 2003 F) MA (2019—2020 Æ) 时 SOHO 卫星 对 26~34 nm 极 紫外 辐射 测 
量 值 和 Sio 指数 


气 密 度 值 (CHAMP 和 GRACE-A 的 大 气 密度 值 由 美国 Colorado 大 学 提供 ，SWARM-C 大 
气 密度 值 由 ESA 任务 数据 中 心 发 布 )。 资 料 时 间 范 围 如 图 1 a) 所 示 ， 这 3 颗 星 资料 总 跨度 


是 2001—2021 年 共 20 a, Mite f 2 次 太阳 活动 峰 年 ，2 KAE. CHAMP 和 GRACE-A 
接近 6 a 的 时 间 重 下，GRACE-A 和 SWARM-C 有 11 月 的 时 间 重 个 ， 可 被 用 于 多 星 资 
LAO EEE. HMA, CHAMP 高 度 在 350~450 km, GRACE-A 高 度 在 460~520 km, 
SWARM-C 高 度 在 480~550 km。 

为 凸显 太阳 辐射 27 d 震荡 对 大 气 密度 的 影响 ， 应 排除 地 磁 显 著 扰动 时 的 资料 ， 选 择 依 
据 是 : 资料 时 刻 对 应 的 日 平均 地 磁 指 数 A, < 20， 且 前 一 天 的 最 大 地 磁 指 数 ap < 30。 此 外 ， 
还 需 识别 在 资料 时 刻 之 前 4 d 内 没有 发 生 3h 指数 ap > 94 的 地 磁 暴 事件 。 如 果 有 ， 则 资料 
很 可 能 处 在 磁 暴 之 后 的 恢复 期 ， 应 被 剔除 。 此 外 ， 偏 心率 引起 轨道 高 度 的 周期 变化 ， 以 及 半 
长 轴 长 期 衰减 均 会 引起 大 气 密度 变化 ， 为 避免 这 些 信 号 掩盖 太阳 辐射 的 影响 ， 需 要 对 不 同 高 
度 资料 做 归 算 处 理 ， 先 计算 研究 区 间 内 的 轨道 高 度 平均 值 ， 再 利用 式 (1) 将 大 气 密度 统一 归 


ao 
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算 到 该 平均 高 度 2 LE: 


= Pm(Z) 
PESE T , (1) 
其 中 ，p(2) 为 高 度 2 处 的 实测 密度 ，py,(2Z) 和 pm(Z) 为 用 模型 计算 的 两 个 高 度 Z TZ 处 
的 大 气 密度 ，p(2) 是 实测 密度 值 归 算 到 Z 处 的 密度 值 。pw 选择 NRLMSISE00 模型 进行 
计算 。 
2.4 ”提取 大 气 密度 与 辐射 短期 变化 的 定量 关系 

主流 模型 一 般 把 大 气 密度 的 变化 区 分 为 随 太 阳 辐 射 的 长 期 变化 和 短期 变化 两 部 分 : 长 
期 变化 用 81 d 平均 射电 辐射 流量 表征 ， 短 期 变化 用 射电 辐射 流量 日 值 减 去 81 d 平均 射电 辐 
射流 量 刻画 。 本 文 按照 相同 的 原理 ， 建 立 大 气 密度 的 简化 回归 方程 式 (BIR (2))， 采 用 极 紫 
外 辐射 指数 S10， 长 期 变化 部 分 用 S. 表征 ， 短 期 变化 用 式 中 第 二 项 (日 值 减 去 滑动 均值 ) 描 
述 ， 本 文 重 点 关注 这 一 项 的 偏 相 关系 数 Vor 和 斜率 Barz?， 式 中 第 三 项 是 地 磁 活 动 对 大 气 密度 


影响。 


— 


p = 01981 + a2(Saaily Sg) ta3kp+a4 . (2) 

对 每 一 天 的 密度 资料 ， 先 根据 星 下 点 地 方 恒星 时 区 分 成 日 照 区 和 阴影 区 两 组 ， 再 对 南北 

b 75° 之 间 的 密度 求 日 平均 值 ， 从 而 研究 明 夜 和 高 度 差 异 ， 在 研究 纬度 差异 时 进一步 按照 低 
EHE (0*~30°)、 中 纬度 (30"~60?)、 高 纬度 (60°~75°) 统计 密度 均值 。 


Mi 
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3 ”大气 密度 27 d 短期 变化 规律 


3.1 ， 高度 与 昼夜 差异 

为 展示 27 d 变化 特征 ， 先 从 大 气 密 度 和 Sio 指数 时 间 序 列 中 扣除 长 期 趋势 项 (81 d TR 
动 均值 )， 剩 余 量 即 为 短期 变化 信号 。 以 大 气 密度 的 时 间 序 列 为 例 ， 提 取 短 期 信号 的 方法 是 : 
W pa 为 大 气 密度 日 均值 ，pr 为 81 d 平均 值 所 描述 的 长 期 变化 ，ps 为 周期 约 为 27 d 的 短 
期 变化 ，t 为 计算 日 期 ， 有 : 


t 十 40 
pr 的 = M palt), ps(t) = palt) — pr(t) . (3) 

t—40 
WAH t. 取 其 前 后 各 40 d (共计 81 d) 计算 均值 ， 随 着 t 的 变化 ， 每 一 步 都 重新 选取 81 d 
窗口 ， 即 窗口 随 着 计算 时 刻 t 滑动 。 由 于 是 事后 研究 ， 辐 射 和 密度 数据 都 可 获取 。 
图 3 a), c) 给 出 CHAMP 和 GRACE-A 在 2004 年 大 气 密度 与 辐射 短期 变化 曲线 ， 两 者 
同步 涨 落 且 幅度 一 致 ，b), d) 给 出 两 者 之 间 相 关系 数 和 斜率， 相关 系数 为 0.897 (CHAMP) 
和 0.847 (GRACE-A)， 呈 现 显著 的 线性 相关 。 和 斜率 的 物理 意义 是 27 d 短期 变化 量 每 变化 一 
个 单位 ， 引 起 大 气 密度 的 变化 量 ， 表 征 了 大 气 密 度 受 辐射 影响 的 程度 ，CHAMP 为 0.0515， 
GRACE-A 为 0.0139， 后 者 是 前 者 的 27%。 由 于 两 颗 星 的 资料 在 同一 时 间 区 间 ， 但 处 于 不 
同 高 度 (2004 年 CHAMP 高 度 均值 390 km, GRACE W 489 km)， 这 意味 着 随 着 高 度 增加 ， 
辐射 对 密度 影响 程度 会 减 小 。 
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3 AIMES Sio 指数 和 日 平均 大 气 密度 之 间 的 关系 (2004 年 ， 日 照 区 ) 


用 更 长 时 间 的 资料 ， 依 托 多 元 回归 方程 ( 见 式 (2)) 统计 得 到 表 1， 给 出 大 气 密度 与 辐射 


27 d 短期 变化 的 偏 相关 系数 Voz 和 回归 系数 Bor (意义 等 同 于 图 3 中 的 斜率 )， 无 论 在 日 


照 区 


或 是 阴影 区 ， 大 气 密度 与 Sio 指数 均 存 在 显著 相关 ( Voz 超过 0.8)。 对 比 2014 年 SWARM-C 
和 GRACE-A 的 资料 ， 同 样 证 明了 高 度 越 高 ， 斜 率 越 小 。 


进一步 对 比 表 1 中 和 白天 和 夜间 半球 的 结果 ， 可 以 看 出 太阳 辐射 对 密度 影响 存在 昼夜 差 


异 : 夜间 密度 与 辐射 的 相关 性 依然 显著 ,但 在 数值 上 Voz 和 Bor 均 比 白天 半球 低 ， 夜 间 的 


Boz 只 有 白天 的 46%~55%， 意 味 着 太阳 辐 冉 短期 震荡 对 夜间 大 气 密 度 的 影响 程度 ， 约 为 白 


天 的 一 半 。 这 主要 是 因为 白天 半球 直接 受到 辐射 影响 产 4 


来 自白 天 半球 能 量 和 温度 的 传播 。 


表 1 中 还 给 出 了 大 气 密度 模型 NRLMSISE00 对 应 的 模拟 结 
结果 始终 低 于 观测 统计 ， 表 明 模 型 低估 了 太阳 辐射 短期 变化 对 密度 的 影响 ， 以 2003— 2004 


E 密 度 涨 落 ， 而 夜间 半球 是 间接 受到 


o (1) 从 B27 看 ， 模拟 


^F. CHAMP 的 日 照 区 资料 为 例 ， 模 拟 斜率 为 0.0221， 观 测 结果 为 0.0543， 模 拟 仅 为 实测 的 


41%; (2) 从 Bor 的 公 夜 比值 看 ， 模 型 与 实测 基本 一 致 ， 


表明 模型 对 昼夜 差异 的 刻画 是 准确 
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表 1 相同 时 期 不 同 高 度 卫星 测量 的 大 气 密度 和 Sio 指数 相关 性 及 模型 模拟 结果 
观测 统计 模拟 结果 
时 间 卫星 资料 高 度 照 区 阴影 区 B27 EK WK B27 
/km | Voy B27 V27 B27 N/D B27 B27 N/D 
CHAMP 400 | 0.924 0.0543 0.890 0.0297 55% | 0.0221 0.0133 60% 
2003—2004 | GRACE-A | 490 | 0.876 0.0147 0.811 0.0068 46% | 0.0068 0.0030 44% 
Diff. 90 27% 23% 31% 23% 
GRACE-A | 430 | 0.884 0.0300 0.841 0.0165 55% | 0.0181 0.0104 57% 
2014 SWARM-C | 480 | 0.878 0.0142 0.878 0.0067 47% | 0.0100 0.0049 49% 
Diff. 50 47% 41% 55% 47% 
注 : 表 中 Bor N/D 表示 用 夜间 斜率 和 日 间 斜 率 求 比值 ，Di 为 两 颗 星 之 间 结 果 的 差异 ， 对 于 高 度 ALT 为 绝对 差 值 ， 
对 于 斜率 B27 为 两 颗 星 之 间 的 比值 ， 比 值 越 小 意味 着 差异 越 大 ， 限 于 NRLMSISE00 模型 的 固有 算法 ， 模 拟 中 采 
Fyo.7 指数 描述 太阳 辐射 水 平 (与 Sio 指数 量 纲 一 致 )。 模 拟 结果 的 相关 系数 Var 始终 超过 0.99， 这 里 不 再 列 出 。 
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AS 


味 着 模型 


LIKE 
3.2 


1 


E1 


13 2011 4 


对 


年 1 月 的 资料 ，5 
可 见 ， 强 太阳 辐射 水 平 


时 ,1 


倍 (白天 半球 )， 意 味 着 


XP 


活动 剧 


烈 时 ，27 d 短 


o 均值 为 142.616。 这 两 段 资料 在 相同 高 度 范围 
局 相关 系数 Vor EK, RR Bo. 是 弱 太 PB 
响 更 显著 。 夜 间 半 球 的 规律 是 


期 


ERI 


Ie 1g 


D 星 做 对 比 研 究 。 表 2 第 1 组， 选择 GRACE-A 1 
KIEL, Sio 均值 为 77.966; 选择 SWARM-C 卫星 2014 年 4 月 至 2015 


， 但 Sio 指数 均值 相差 65, 


H DTM2000 模型 ; 


(3) 对 比 不 同 高 度 的 Bor, Œ 2014 年 模拟 的 比值 (5596, 47%) 普遍 比 观测 (47%, 41%) 
定 程度 上 低估 了 高 度 增加 带 来 的 斜率 衰减 。 我 们 还 月 
， 结 果 与 NRLMSISE00 模型 一 致 。 
不 同 辐射 水 平时 的 差异 
选择 高 度 范围 相同 的 两 颗 ] 


行 


HE 2010 年 


IA 
HZ 


$ 


SR INS EA 1.760 


致 的 。 对 于 第 2 4H, SWARM-C RX 2016 年 4 月 一 2017 4E 3 月 的 资料 ，GRACE-A 取 为 


2012 年 5 月 一 
Hj. 辐射 水 平 越 高 ，27 d 短期 震荡 的 影响 越 


2013 年 7 


大 。 


表 2 给 出 NRLMSISE00 模型 的 模拟 结果 作为 参 
以 第 一 组 GRACE-A 的 日 间 资 料 为 例 ， 模 拟 斜率 


于 观测 统计 ， 


(0.0084) 的 49%, SH 1 中 结果 一 致 ， 进 
(2) 从 斜率 比值 看 ， 模 拟 结 果 
表 1 和 表 2 只 是 针对 个 别 
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时 间 


同 ”， 
(按时 i 
系 。 
是 夜 
度 范 


SB 


[uni 
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接近 )， 


以 及 “高 度 相同 、 委 
司 连 续 分 为 16 组 ) 统计 
图 4 可 知 : (1) 斜率 随 高 度 增 大 而 衰减 ， 
my 1.85 fi: (2) E 4 c), d) KEET GRACE-A 和 SWARM-C 的 结 


42H 


F 


到 斜率 Boz, 


斜率 随和 


差异 


开始 显现 ， 同 时 ， 这 两 


昌 射 增强 而 增 大 ， 呈 了 


EAE" PRRHBUE FÉ 
并 绘制 Bor 与 资料 所 处 的 高 度 、 辐 
4a) Rib), E 


步 而 


的 资料 ， 两 段 资料 高 度 均值 相同 ， 辐 射 


照 。(1) AR 


RE 


均 


Lt 


相差 约 50。 同 样 证 


ES Boz 看 ， 模拟 斜率 小 
(0.0041) 仅 为 观测 结果 


有 认 模 型 低估 了 辐 贡 


短期 变化 对 密度 的 影响 ， 
著 高 于 观测 统计 ， 说 明 模型 高 估 了 由 辐射 水 平 引起 的 差异 。 


EXE: 


EH 


因素 的 作用 趋势 在 


X ( 见 图 


Ko HAMIL, AB 


区 间 内 的 资料 作对 比 ， 分 别 揭示 了 “辐射 相 


同 、 高 度 不 


辐射 影响 差异 。 我 们 还 对 3 颗 星 全 部 资料 


FI 
FH 


fi 


射 水 平 的 关 


(B 


射 对 大 气 密度 的 


Bel 


JAR ERARIK FW, Ee rae REA As ES, Dm 


" 


天 半球 均 比 夜间 


[ 


接近 时 ， 辐 射 水 3 
球 显著 。 


间 衰 减速 度 
为 两 颗 星 高 
影响 


区 
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表 2 相同 时 期 不 同 高 度 卫星 测量 的 大 气 密度 和 Sio 指数 相关 性 及 模型 模拟 结果 


观察 统计 模拟 结果 

a E 日 照 区 阴影 区 日 照 区 ”阴影 区 

按 高 度 分 组 卫星 资料 S10 均值 Vay By Vo By Bor By 
Group 1 GRACE-A 77.966 0.681 0.0084 0.431 0.0025 | 0.0041 0.0016 
Z= ATI km SWARM-C 142.616 | 0.899 0.0148 0.868 0.0073 | 0.0102 0.0050 

Diff. 64.650 1.76 2.87 2.49 3.13 
Group 2 SWARM-C 77.978 0.777 0.0128 0.732 0.0057 | 0.0056 0.0023 
GRACE-A 128.428 | 0.935 0.0162 0.867 0.0074 | 0.0112 0.0057 

Z = 458 km 


Diff. 50.450 1.26 1.29 2.00 2.48 
YE: KF Diff. 为 两 颗 星 /两 组 资料 之 间 的 差异 ， 对 于 Sio 均值 指数 为 绝对 差 值 ， 对 于 B27 为 比值 (中 / 强 辐 
射 情形 除 以 弱 辐射 情形 )， 比 值 越 大 意味 着 差异 越 大 。 模 拟 结果 的 相关 系数 Var 始终 超过 0.99， 这 里 不 再 列 出 。 
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4 3 颗 卫 星 16 组 资料 分 别 统计 得 到 斜率 Bor 随 高 度 和 辐射 水 平 的 变化 规律 


3 期 


3.3 


K 3 给 出 4 组 资料 按照 
发 现 : 


大 人 


纬度 差异 


天 和 夜间 半球 分 别 统计 的 3 个 纬度 范围 
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409 


的 Voz 和 B27。 对 比 可 


在 日 照 区 ， 和 斜率 随 着 纬度 增 大 而 减 小 ， 低 纬度 区 域 是 斜率 最 大 值 ， 高 纬度 区 域 斜 率 是 
低 纬度 的 80%~85%; 在 阴影 区 情形 相反 ， 和 斜率 随 着 纬度 增 大 而 增 大 ， 高 纬度 区 域 是 斜率 最 
， 高 纬度 是 低 纬 度 的 1.1~1.3 倍 。 
表 3 3 颗 卫 星 不 同 纬度 资料 统计 得 到 的 Voz 和 Boz 
a " 照 区 阴影 区 

卫星 资料 和 时 间 范 转 纬度 V Bu V Bu 

高 纬 | 0.890 0.05478 | 0.867 0.036 33 

CHAMP 中 纬 | 0.894 0.06414 | 0.861 — 0.03407 

2003—2005 KZ | 0.914 0.06518 | 0.870 0.03339 

高 / 低 续 0.840 1.088 

高 纬 | 0.838 0.01082 | 0.783 0.007 25 
GRACE-A 中 纬 | 0.849 0.01248 | 0.782 0.00557 

2003— 2005 KA | 0.852 0.01305 | 0.799 0.00545 

高 / 低 续 0.829 1.330 

高 纬 | 0.786 0.01422 | 0.759 0.00888 

GRACE 中 纬 | 0.807 0.01685 | 0.730 0.00795 

2011 一 2013 Kai | 0.828 0.01664 | 0.711 0.00713 

高 / 低 续 0.855 1.245 
高 纬 | 0.794 0.01459 | 0.757 0.00915 
SWARM-C 中 纬 | 0.801 0.01704 | 0.772 0.00802 
2014—2016 Kai | 0.811 0.01723 | 0.780 0.00726 

高 / 低 续 0.847 1.260 


4 大 气 密度 对 极 紫 外 辐射 的 啊 应 时 间 


北纬 
射 测 


义 上 


构建 回归 方程 


13d 


选择 1 h 43 


75° 之 间 的 平均 大 气 密度 ， 得 到 的 密度 时 序 资料 分 辨 率 为 3 h 
量 值 和 其 他 影响 因素 (地 磁 、 地 方 时 等 ) 
(72~24 h)， 观 察 回归 方程 相关 系数 最 大 (回归 效果 最 优 ) 所 对 应 的 延迟 时 间 ， 即 为 统计 意 


SN 度 对 太阳 极 紫 外 辐射 的 响应 时 间 。 


利用 1h 分 辨 率 的 SOH 


Vet 


的 多 元 回归 方程 ， 通 过 


PRAI 2634 nm 极 紫 外 测量 值 ， 再 计算 3 个 连续 上 升 (或 下 降 ) 轨道 南 


过 改变 延迟 


HE KA EE 55 58 


SHY Ta FI TL 


O 卫星 SEM 测量 值 ， 与 3 h 分 辩 率 的 大 气 平均 密度 (临近 3 
， 见 式 (4 和 )。 为 凸显 大 气 的 密度 短期 变化 ， 以 


为 一 个 计算 窗口 ， 因 变 


量 选择 为 大 气 密度 


(1) 极 紫外 辐射 的 短期 变 
(2) 地 磁 指 数 Kp 的 短 其 


化 量 AFguy; 
HEE AKp; 


(3) 地 方 时 的 变化 量 ALT; 


的 变化 值 Ap。 回 归 方 程 的 自 变 


量 考虑 6 类 : 
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(4) 地 方 时 LT 的 传 里 时 函数 形式 ; 
(5) 计算 窗口 [ti ti +13 d] 极 紫外 辐射 均值 Exvv: 


(6) 计算 窗口 内 平均 轨道 高 度 Z. 


从 各 


因素 的 作用 原理 上 看 ，(1) 一 (3) 类 属于 在 同一 个 27 d 周期 内 ， 大 气 密度 变化 受到 


的 直接 影响 ， 由 于 3 颗 卫 星 的 极 轨道 特点 ， 在 计算 窗口 内 地 方 时 变化 量 ALT, 5 最 
大 值 小 于 2 h， 所 以 第 三 项 ALT 与 Ap 的 关系 可 近似 为 线性 。(4) 一 (6) 类 属于 不 同 的 27 d 


的 周期 之 间 影 响 大 气 密度 变化 的 
表示 ， 与 2 的 关系 可 以 近似 为 线性 关系 。 


辐射 相对 变化 


具体 建立 样本 的 方法 如 示意 图 
选取 的 延迟 时 间 ，5 为 计算 窗口 内 的 平均 大 气 密度 ， 有 : 


4) — pl) fp 
Eruv(tij + tag) — Eguv(ti + tiag)|/EEuv 
tij) — Kp(ti) 


因素 ， 密 度 随 地 方 时 LT 的 周 日 变化 可 以 用 2 阶 传 里 叶 函 数 


Ap = b; AE guy + bo, AK p + bz ALT + ba cos(LT) 十 
bs sin(LT) + bg cos(2LT) + b; sin(2LT) + bsEgov + 09Z.. (4) 


5, ilti; — t; +j (j =7~ 13), tag 为 计算 Bpuv ANA 


Lnd 


图 5 


时 间 t/d 


计算 窗口 内 样本 选择 示意 图 


ae 


3 期 EER 
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， 等 ; 热 层 大 气 密度 对 太阳 辐射 27 d 短期 震荡 的 响应 规律 


按照 时 间 序 列 ， 将 去 之 后 的 13 d 作为 
对 应 求 差 。 为 凸显 极 紫 外 辐 


射 变化 A Eguv(i,;) J) 的 影响 , 


—^ BUE, H t; 时 刻 的 资料 分 别 与 t; 时 刻 资料 
还 应 尽量 降低 地 磁 扰动 ACKC p; ; 对 密 


度 的 影响 ， 所 以 剔除 Kp(t;) > 4, Kp(ti;) > 4, AKpi;; > 2 的 资料 。 计 算 Kp 时 选择 资料 时 
口 ， 得 到 式 (4) 


刻 前 3 h 和 6 bh 的 地 磁 指 数 加 权 平 均 。 通 过 


中 参数 估计 所 需 的 样本 。 
在 一 


SEA ts:， 可 以 得 到 大 量 的 计算 窗 


度 延 迟 响 应 。 


EIKES 
a) PAS 


变化 量 的 1 
的 是 辐射 上 升 阶段 的 情况 ， 


为 保证 资料 的 充分 性 ， 我 们 出 


为 阴影 


ba 3 颗 卫 星 的 资料 对 回 
相关 系数 如 


区 大 气 并 不 直接 受到 太阳 辐射 作 月 


—18 h, 


Bi 18 h 


6 A Hanae 况 ， 此 时 1 


峰值 ， 


0 为 极 紫 外 辐 
吕 射 的 响应 时 


归 方程 进行 求解 ， 
图 6 a) 和 c) 所 示 ， 横 坐标 是 选择 不 同 的 延迟 时 间 : 
偏 相 关系 数 最 大 值 出 现在 
的 极 紫外 辐射 描述 大 气 密度 短期 震荡 ， 回 归 效 果 最 好 ， 最 大 相关 系数 为 0.897; 
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E 


大 气 密度 变化 


图 6 


意味 着 用 密度 时 刻 之 
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一 36 h 都 相对 接 


峰值 附近 比较 平坦 。 图 6 b) 和 d) 是 选取 上 升 段 结果 -18 hn， 下 降 段 结果 -30 hp， 计算 
pes 度 变 化 量 与 极 紫 外 变化 量 的 散 点 图 ， 可 以 看 出 清晰 的 线性 规律 。 
O90 g ama en L2 B EIDE 
P4 N. | 4 =18 h 
0.85} f,,7 18 h n 
d P4 R „=0.897 N 0.8} 
m 0.80} , X E 
E 0.75+ J \ 4 0.A- ; 
0.70} / i? E 
0.65 =72 -48 -24 0 24 0.00 0.10 0.20 0.30 0.40 
t, /h AE Buy 
a) b) 
0.85} Euv s FEB pan 0.0} Fan rior 
ger 
Zo qu 30h 
wa 0.754 7 R „=0.852 \ -0.4 
4K Le N 4 
E 0.65| WE -0.8 
0.55L LU SER E HERE NE 
=72 —48 -24 0 24 -0.35 -0.25 -0.15 -0.05 
ts/ AE guy 
c) d) 
6 ”大 气 密度 对 极 紫 外 辐射 的 响应 延迟 时 间 (区 分 辐射 上 升 段 、 下 降 段 ) 


412 
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归 效 果 略 低 于 上 升 段 ， 且 延迟 时 间 晚 于 上 升 段 ， 是 
极 紫外 辐射 输入 减弱 ， 作 用 下 降 ， 而 凸显 出 地 磁 活 动 的 影响 相 
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归 效 果 变 差 ， 峰 值 附近 变 得 
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在 物理 大 气 模型 TIEGCM 模拟 中 将 Kp 设置 为 常数 值 0， 单 独 分析 极 紫外 辐射 对 大 气 密度 
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5 小结 与 讨论 


本 文选 择 CHAMP、GRACE-A £l SWARM-C 共 3 颗 极 轨 卫 星 获 
并 从 中 提取 大 气 密度 的 27 a 短期 变化 信号 ， 然 后 与 太阳 极 紫 外 辐射 27 d 短期 震荡 时 间 序 列 
进行 相关 性 统计 。 相 比 过 去 的 研究 ， 本 文 使 用 的 大 气 密度 资料 具有 高 时 空 分 辨 率 和 全 球 纬度 
的 太阳 辐射 资料 (S10 指数 和 SEM 测量 值 ) 是 直接 驱动 大 气 密度 变化 的 极 
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因素 ， 辐 射 水 平 次 之 。 将 观测 结果 与 
辐射 短期 变化 对 大 气 密度 影响 的 
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响应 (延迟 ) 时 间 ， 约 为 18~30 
意义 上 确定 延迟 时 间 为 18 h。 


将 观测 特征 与 模型 模拟 结果 对 比 ， 表 明 现 有 模型 的 建 模 尚 有 不 足 。 后 


和 DTM 模型 框架 作 建 模 试验 ( 


(3) 纬度 差异 ， 白 天 半 球 低 纪 


度 区 域 受 到 辐 
形 相反 ， 夜 间 高 纬度 地 区 是 


响 程 度 逐 步 降 低 。 考 虑 到 能 量 在 传播 过 程 中 会 不 断 耗 散 ， 以 上 纬度 差异 可 能 预示 着 白天 至 夜 
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能 量 传播 的 路 径 是 : 先 从 白天 中 低 纬度 向 高 纬度 传播 ， 再 到 达 夜 间 高 纬度 区 域 ， 随 
全 球 大 气 的 整体 涨 落 。 

(4) 密度 对 辐射 的 响应 时 间 : 以 白天 半球 为 例 ， 厂 
h。 考 虑 到 处 于 辐射 上 升 段 的 资料 信 噪 比 更 强 ， 
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续 拟 采用 Jacchia 
， 从 外 层 温 度 与 辐射 的 线性 关系 式 入 手 ， 对 于 短 


期 变化 项 (日 值 与 81 d 均值 2 
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测 中 的 特征 ， 积 累 数值 试验 可 以 优化 建 模 。 


重新 计算 偏 导数 和 模型 系数 ， 观 察 能 否 有 效 提取 观 
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Abstract: In the mission of orbit determination and prediction for space object and debris in 
low earth orbit, accurate mass densities and its short-term variation of thermospheric neutral 
gas are required. Solar radiation’s short-term variation with 27-d period is an important 
space environment factor because it drives the global-scale fluctuation of thermospheric mass 


densities synchronously. This paper uses the neutral mass densities data derived from the 
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non-gravitational acceleration measurements by three polar-orbit satellites, such as CHAMP, 
GRACE-A and SWARM-C. The short-term variation of these densities was extracted and 
compared with the solar extreme ultraviolet (EUV) radiation, indicated as S10 index, to 
study the thermosphere response to radiation. Their differences at various height, radiation 
level, dayside/nightside hemisphere, and latitude were discussed. The multivariate regression 
analysis method was used to calculate the correlation coefficient and slope between the 
density and solar radiation. The results were as follows. (1) solar variation’ effect on the 
density is negative correlation with altitude (major factor) and positive with solar activity 
level (minor factor). (2) the effect in dayside hemisphere is about 1 time stronger than that 
in the nightside. (3) the effect decreased as the rising of latitude in the dayside, while in the 
nightside the opposite is the case. It is implied that the direction of energy and temperature 
conduction in the dayside is from low equator region to high latitude, and in the nightside 
from high to low latitude. By comparing these observation results with the simulation 
used NRLMSISE00 model, it is indicated model underestimate the effect of solar radiation’s 
short-term variation on densities for about 50%~60%, and its difference at various altitude, 
while overestimate the difference caused by average level of solar radiation. The response 
lag time of densities to solar radiation was studied by using EUV measurement with 1-h 
temporal resolution from SOHO project, indicating the lag time was about 18 h (0.75 d). 
These conclusions showed the complicate effect of solar radiation’s short-term variation and 


were helpful to improving the modeling atmosphere density. 
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